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On  the  basis  of  the  theory  of  a-stability,  developed  in 

a  previous  paper,  a  classification  of  high-shear  o-stable 

diffuse  pinch  configurations  is  obtained  by  numerical  solution 

of  the  Euler  equation  corresponding  to  the  modified  energy 

principle.   Toroidal  effects  are  introduced  by  putting  a  lower 

bound  on  the  magnitude  of  the  poloidal  field  in  order  to 

assure  toroidal  equilibrium  and  by  deriving  and  employing  the 

proper  generalization  of  the  Kruskal-Shafranov  limit  for  m  =  l 

unstable  configurations.   The  classification  of  high-shear 

a-stable  configurations  consists  of  five  different  types,  viz. 

Tokamak,  Bg-reversed  pinch,  screw  pinch,  combined  pinch, 

B  -reversed  pinch.   Among  the  significant  new  results  are: 

z 
complete  stability  of  some  of  the  Tokamak  configurations  for  p 

up  to  hio,    a  stable  high-P  [-^^o)   Bg-reversed  pinch  configuration 

that  has  not  been  investigated  before,  a  stable  diffuse  screw 

pinch  with  p  =  30^  when  currents  are  allowed  to  flow  at  the 

wall  and  with  p  =  8^  with  vanishing  current  at  the  wall,  the 

possibility  of  containing  diffuse  high-p  m  =  l  stable  plasmas 

without  B  -field  reversal  or  pressure  minimum  on  axis.   Here, 

a-stability  refers  to  times  needed  for  fusion. 


1.    Introduction. 

In  a  precedir 

introduced  the  concept  of  o-stability,  observing  that  for 
thermonuclear  confinement  of  plasma  it  is  not  necessary  to 
demand  absolute  stability  but  that  it  is  sufficient  that 
perturbations  do  not  grow  faster  than  e   when  1/a  corre- 
sponds to  some  characteristic  time  needed  for  fusion.   In 
this  paper  the  stability  of  the  diffuse  linear  pinch  is 
reinvestigated  on  the  basis  of  this  enlarged  stability 
concept,  where  the  aim  is  to  derive  a,  as  complete  as 
possible,  classification  of  a-stable  equilibria.   The  sig- 
nificance of  this  approach  rests  upon  the  fact  that  in 
ideal  magnetohydrodynamics  many  unstable  modes  occur  which 
have  vanishingly  small  growth  rates,  so  that  the  class  of 
o-stable  equilibria  which  are  acceptable  in  fusion  devices 
is  much  wider  than  the  class  of  absolutely  stable  equilibria. 
In  demanding  absolute  stability  one,  therefore, dismisses 
a  number  of  interesting  experimental  possibilities  which 
could  have  advantages  other  than  following  from  ideal  MHD . 

In  the  usual  stability  analysis  of  the  diffuse  linear 
pinch  a  local  criterion  (Suydam's    criterion)  is  fairly 
easily  obtained,  but  it  is  well-known  that  local  stability 
does  not  guarantee  global  stability.   In  order  to  establish 
global  stability  one  has  to  solve  a  differential  equation, 
where  all  the  details  of  the  equilibrium  distribution  of 
magnetic  field  components  and  pressure  enter.   Stability 


then  often  turns  out  to  be  such  an  extremely  sensitive  function 
of  the  equilibrium  profiles  that  the  relevance  for  experimental 
situations  becomes  very  doubtful.   For  this  reason,  a  o- 
stability  analysis  is  much  more  satisfactory,  because  many  of 
the  subtleties  of  ideal  MHD  connected  with  the  singular 
behavior  of  the  equations  at  marginal  stability  disappear. 
Consequently,  it  becomes  possible  to  construct  bands  of  a-stable 
equilibrium  profiles,  so  that  a  certain  amount  of  imprecision 
is  allowed  for  a  profile  in  order  to  be  a-stable.   Such  an 
approach  seems  to  be  a  proper  theoretical  reflection  of  the 
fact  that  experimentally  profiles  cannot  be  produced  with  any 
desired  amount  of  precision. 

The  present  paper  is  entirely  devoted  to  a  numerical 
investigation  of  a-stable  diffuse  linear  pinch  configurations. 
The  analysis  is  based  on  the  oscillation  theorem  for  a- stabil- 
ity derived  in  I,  Sec.  3a.   This  theorem  is  analogous  to 
Newcomb's-^  "^  theorems  pertaining  to  marginal  stability,  but  it 
has  the  important  advantage  that  no  singularities  are  encoun- 
tered so  that  a  faster  numerical  calculation  is  possible.   Thus, 
computer  time  is  saved  for  the  purpose  of  scanning  many  equilib- 
rium profiles.   The  only  restriction  on  the  form  of  the  chosen 
equilibrium  profiles  has  been  that  it  should  satisfy  Suydam's 
criterion.   The  implication  of  this  restriction  has  been  dis- 
cussed at  length  in  I,  Sec.  4.   For  the  present  paper,  it 
implies  that  we  are  restricting  the  analysis  to  intermediate- 
and  high-shear  configurations  where  violation  of  Suydam's 
criterion  in  most  cases  leads  to  rapidly  growing  nonlocalized 


m =  1  modes  having  large  growth  rates.   Therefore,  Suydam's 
criterion  provides  a  convenient  way  of  preselecting  equilib- 
rium profiles  for  the  intermediate-  and  high-shear  cases. 

It  should  be  clear  from  the  preceding  discussion  that 
a-stable  diffuse  pinch  configurations  are  to  be  found  in  two 
broad  categories,  viz.  high-shear  configurations  and  low- shear 
configurations.   The  first  class  is  treated  in  this  paper, 
where  the  returning  feature  will  be  improved  stability  by 
increasing  the  shear.   The  existence  of  the  second  class  of 

low-shear  a-stable  configurations  has  been  indicated  by 

h.) 
Grad.  '      It  will  be  the  subject  of  future  numerical  investiga- 
tions.  The  reason  for  a  separate  treatment  is  that  the  rules 
for  obtaining  and  optimizing  the  two  types  of  configurations 
are  completely  different.   For  the  high-shear  configurations 
Suydam's  criterion  acts  as  a  necessary  condition  for  a-stabil- 
ity  (except  in  a  few  cases  to  be  found  in  the  text).   Toroidal 
periodicity  is  taken  into  account  only  for  m=l,  the  unstable 
regions  for  higher  m  normally  being  broad  and  sometimes  over- 
lapping.  For  the  low- shear  configurations,  on  the  other  hand, 
Suydam's  criterion  is  not  relevant  and  the  higher  m  modes  are 
stabilized  by  toroidal  periodicity.   In  this  case,  a-stability 
provides  the  connection  between  favorable  closed-line  stability 
results  and  the  possibility  of  low-shear  a-stable  systems. 

Most  of  the  older  investigations  on  the  linear  pinch  have 
been  restricted  to  sharp-boundary  and  constant-pitch  models. 
With  the  exception  of  Ref .  6,  these  studies  have  only  more 
recently  been  extended  to  include  numerical  calculation  of 
growth  rates  of  instabilities  for  completely  diffuse 
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models'^^'®  ^'^  .   One  of  the  merits  of  the  numerical  approach 
is  that  one  can  easily  check  whether  the  results  obtained  from 
special  models  have  a  more  general  validity.   Although  a  number 
of  results  do  stand  this  test,  there  is  also  a  fair  amount  of 
deep-rooted  misconceptions  based  on  the  results  obtained  from 
special  models,  e.g.:   that  the  stability  of  the  ra=l  mode  does 
not  very  much  depend  on  the  details  of  the  current  distribution, 
that  violation  of  Suydam's  criterion  leads  only  to  localized 
modes,  that  the  ©-pinch  is  a  singularly  isolated  stable 
equilibrium,  that  sharp-boundary  pinches  are  more  stable  than 
diffuse  ones,  or  vice  versa,  that  one  can  distinguish  between 
kinks  and  interchanges,  that  the  worst  perturbations  are 
incompressible,  that  incorapressibility  is  a  good  approximation 
if  the  longitudinal  field  is  large,  that  the  generalization  of 
the  Kruskal-Shafranov  limit  for  a  diffuse  pinch  can  be  expressed 
in  terms  of  a  local  criterion  for  q,  etc.,  etc.   It  is  an  easy 
matter  to  produce  for  each  of  these  statements  a  profile  which 
contradicts  it. 

The  main  problem  in  abstracting  general  information  on 
the  question  about  what  is  good  and  what  is  bad  for  stability 
is  the  fact  that  stabilizing  and  destabilizing  factors  cannot 
be  separated.   E.g«,  a  current  parallel  to  the  magnetic  field 
drives  instabilities  but  also  creates  stabilizing  shear,  so 
that  the  stability  just  depends  on  the  details  of  the  current 
distribution.   Our  approach  to  this  problem  here  is  to  generate 
as  many  different  a-stable  profiles  as  possible,  in  this  way 
abstracting  information  on  this  question. 


This  paper  is  arranged  as  follows.   In  the  next  section 
the  basic  equations  are  given  and  the  ninnerical  procedure  is 
outlined.   In  Sec.  3  simple  conditions  for  toroidal  equilibrium 
and  stability  are  introduced  as  additional  constraints  in  the 
theory.   In  Sec.  k,    which  constitutes  the  core  of  this  paper, 
the  obtained  o-stable  equilibrium  profiles  are  shown.   These 
equilibria  include  Tokamak,  Bg-reversed  pinch,  screw-pinch, 
combined  pinch,  B  -reversed  pinch,  and  different  hybrid  forms 
of  these  configurations.   Roughly  two  categories  are  obtained 
depending  on  whether  or  not  toroidal  periodicity  has  to  be  taken 
into  account.   In  order  to  obtain  a  complete  survey,  we  have  not 
hesitated  to  include  known  results  as  well.   In  Sec.  5  some 
very  unstable  but  realistic  configurations  are  also  shown, 
serving  as  a  warning  against  premature  optimism. 


2.    Procedure. 

We  investigate  numerically  the  a-stability  of  diffuse 
linear  pinch  configurations.   The  equilibria  are  described 
by  the  three  profiles  p(r),  B  (r) ,  and  B  (r) ,  which  are 
connected  through  the  equilibrium  equation 

(p  +  ^  B^)'+   B^/r  =  0.  (1) 

Fixing  the  value  of  B,  where 

B  =  E(0) ,  (2) 

p{0)  +  B^(0)/2 
this  leaves  two  of  the  three  profiles  to  be  chosen  arbitrarily, 
e.g.,  p  and  v   =  B  /rB  .   In  order  to  reduce  the  number  of 

possibilities  further  we  restrict  the  analysis  to  those 
equilibria  for  which  the  pressure  is  a  monotonically  decreas- 
ing function  of  r  (p'<0)  and  which  are  Suydam-stable.   The 
first  restriction  is  a  natural  one  in  view  of  the  fact  that 
we  are  interested  in  confined  plasmas,  whereas  the  second 
restriction  was  discussed  at  length  in  I,  Sec.  4.   From 
Suydam's  criterion 

p.  ,  Ir  b2(J')'  >  0  (3) 

we  now  have  that  |y'/p^l  t'  0  everywhere  if  p'  <  0,  so  that  u 
also  has  to  be  a  monotonic  function  of  r. 

This  leads  to  the  basic  possibilities  depicted  in 
Fig.  1.   The  function  y  can  be  a  positive  decreasing  function 
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of  r  (Fig.  lA) ,  it  can  change  sign  when  B  passes  through 
zero  (Fig.  IB) ,  it  can  be  a  positive  increasing  function 
(Fig.  1C,D) ,  or  it  can  go  to  infinity  and  change  sign  when 
B  passes  through  zero  (Fig.  IE) .   The  cases  C  and  D  will 
be  distinguished  later  on  the  basis  of  toroidal  considera- 
tions.  For  convenience,  we  have  attached  familiar  names 
to  these  configurations  which  does  not  necessarily  imply 
that  the  experimental  confinement  systems  known  by  these 
names  do  have  precisely  the  depicted  profile.   We  shall 
show  that  in  each  of  the  indicated  group  of  profiles  of 
Fig.  1  a-stable  configurations  exist  with  acceptable 
values  of  a  and  e  as  regards  to  fusion.   Many  of  the  a-stable 
configurations  turn  out  to  be  smoothly  transformable  into  one 
another,  so  that  a  lot  of  hybrid  forms  can  be  obtained.   Only 
the  most  basic  types  will  be  presented  here. 

Different  equilibrium  profiles  are  generated  and  tested 
with  respect  to  o-stability.   The  numerical  procedure  is 
based  on  the  theorem  for  a-stability  derived  in  I,  Sec.  3a. 
For  a  certain  range  of  values  of  m  and  k  one  has  to  integrate 
the  a-Euler  equation 


.     [(pg^+F^)   ^po^(Yp+B^)    +   YpF^?l  d_ 
^  L  _  r  dr 


D 

2,(  ^2„2    „„2„2    ^    _^2 

(4) 


D 

2, 


& 


(     kB,      mB^  pa^(Yp+B2)    +YPf2)] 
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where  F   =  mBg/r   +   kB^,  (5) 

24  2222  2  222  2 

D   =    p^a      +    pa    (m''/r''+k'^)  (yp+b'')    +    (m'^/r^+k"')  ypF'' ,  (6) 

subject  to  the  boundary  condition 

rSj.  =  0   at  r  =  0.  (7) 

If  the  obtained  solution  has  no  zero  in  {0,a)  the  pinch  is 
a-stable  for  that  value  of  m  and  k,  if  there  is  a  zero  the 
pinch  is  o-unstable,  whereas  a  zero  at  r  =  a  implies  that 
there  exists  an  eigenvalue  u  =  -a    .      It  is  clear  that  the 
numerical  procedure  can  easily  be  extended  to  calculate 
growth  rates  if  that  is  desirable.   The  range  of  values  of 
m  and  k  to  be  tested  for  a-stability  can  be  reduced  con- 
siderably by  realizing  that  only  the  lower-m  modes  can  be 
unstable  if  Suydam's  criterion  is  satisfied  (I, Sec. 4), 
whereas  the  width  of  the  unstable  k-region  usually  is  only 
of  the  order  | ym | .   This  facilitates  a  rapid  test  for  a-stability. 

The  values  of  o  to  be  chosen  depend  on  the  purpose  of  the 
a-stability  analysis.   For  a  specific  experiment,  a  relevant 

choice  would  be  o  =  /t ,     ,  where  t ,     is  the  available 
decay         decay 

time  as  determined  by  the  decay  of  the  external  or  plasma  currents. 
From  a  a-stability  analysis  one  could  then  find  out  whether  the 
decay  of  plasma  currents  could  be  due  to  MHD  instabilities  or 
should  be  sought  elsewhere.   Here,  we  choose  a  to  correspond  to 
times  needed  for  fusion  according  to  Lawson's  criterion,  so 
that  we  obtain  a  classification  of  diffuse  linear  pinch  con- 
figurations which  are  stable  on  that  time-scale.   The  parameter 


a   enters  the  analysis  only  through  the  dimensionless  quantity 

r  ^  =  i^^£4^  ,  (8) 


which  has  to  be  fixed  by  a  reasonable  choice  of  the  values 

of  p,  a,  and  B  .   For  a  hydrogen  plasma  the  connection 
between  the  e-folding  time  t  and  the  value  of  ^  is  given  by 

the  relation 

1/2 
T  =  4.6  X  10~^^  ^         ^  "^  [sec],  (9) 

o 

where  a  is  given  in  cm,  n  in  particles  per  cm  ,  B   in  Gauss. 
An  example  is  given  in  Table  I  for  a  plasma  with  a  density 
somewhere  in  between  usual  Tokamak  and  pinch  densities.   In 
order  to  distinguish  between  a-stability,  where  growth  at 


0  of  unstable  normal 
modes  e  '^   by  the  dimensionless   quantity  -    Q      =   ^   as 
given  by  Eq.  (8) . 

a-stable  plasmas  of  thermonuclear  interest  are  obtained 
for  values  of  JZ   ranging  from  10~   to  10~   .   It  should  be 
noticed  that  the  precise  value  of  JZ   does  not  enter  the 
analysis  in  a  very  critical  manner.   The  important  feature 
is  just  that  the  mentioned  range  of  values  of  Yl        is  far 


instabilities  destroy  the  plasma  immediately  whereas  it  is 
also  far  away  from  ^   =  0,  where  the  singular  nature  of 
the  marginal  stability  analysis  causes  the  stability  question 
to  be  a  very  critical  one,  depending  sensitively  on  the 
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equilibrium  profiles.   In  the  a-stability  analysis,  on  the 
other  hand,  stable  profiles  occur  in  bands,  so  that  important 
experimental  questions  can  be  answered,  such  as  "How  accurate 
should  a  given  profile  be  produced  in  order  to  be  stable 
during  a  given  time  interval?".   In  Section  5  it  will  be 
shown  that  this  band  is  very  narrow.   However  there  are  an 
infinite  number  of  parallel  a-stable  bands,  hereafter  called 
a-stable  striae  of  equilibrium  profiles,  covering  continu- 
ously the  space  between  two  a-stable  profiles  of  a  same  class. 
There  are  also  a  number  of  a-stable  striae  which  cross  other 
a-stable  striae,  thus  forming  a  mesh  of  o-stable  striae.  An 
example  of  such  a-stable  striae  of  equilibrium  profiles  is 
given  in  Fig.  2,  where  the  pressure  profile  is  kept  fixed  and 
the  B  -profile  is  changed  continuously  from  a  stable  one  with- 
out a  field-reversal  to  a  stable  one  with  a  field-reversal. 
In  between  these  two  cases  a  whole  variety  of  a-stable  fields 
occur. 

From  Eq.  (4)  it  is  evident  that  the  density  profile 
p  (r)  also  enters  the  a-stability  analysis  because  a  only 
appears  through  the  combination  pa   in  the  equation.   We 
are  entirely  free  to  choose  this  profile,  or  the  temperature 
profile,  because  there  is  no  equilibrium  equation  fixing  it. 
However,  the  results  are  very  insensitive  to  the  choice 
of  this  function  (at  least,  as  long  as  it  does  not  become 
very  small) ,  so  that  there  is  little  point  in  considering 
p  (r)  as  an  additional  profile  in  the  classification.   For 


convenience,  we  have  chosen  a  constant  density  throughout 
the  whole  region  (0,a),  which  corresponds  to  a  temperature 
profile   having  the  same  radial  dependence  as  the  pressure 
profile.   The  advantage  of  this  choice  is  that  one  can  easily 
renormalize  the  value  of  YL       if  °^^   wishes  to  take  the 
influence  of  lower  density  regions  into  account.   Moreover, 
we  thus  have  chosen  to  represent  the  outer  regions  of  a 
pinch  as  a  pressureless  plasma  rather  than  a  vacuum, in  this 
way  avoiding  difficulties  associated  with  the  limit  p  -^  0. 
The  latter  difficulties  stem  from  the  fact  that  in  this 
limit  W  -»■  W,  so  that  the  singular  marginal  stability  analysis 
is  recovered.   For  a  real  vacuum  one  can  avoid  singularities 
by  transforming  to  the  variable  Q   (the  radial  component  of 
the  magnetic  field  perturbation) ,  but  in  reality  the  outer 
region  of  a  pinch  is  seldom  properly  approximated  as  a  vacuum 
and  for  some  configurations  this  region  even  carries  large 
equilibrium  currents.   In  this  case,  the  transformation  to 
the  variable  Q   still  makes  sense  if  the  marginal  stability 
with  respect  to  tearing  modes  is  considered,  but  an  investi- 
gation of  tj-stability  with  respect  to  these  resistive  insta- 
bilities would  require  the  study  of  the  full  resistive  MHD 
equations  which  goes  well  beyond  the  scope  of  present  analysis. 


3.    Toroidal  effects. 

A  naive  scheme  for  the  maximization  of  6  with  respect 
to  the  stability  of  a  diffuse  linear  pinch  automatically 
leads  to  the  e-pinch  (B  =  0) .   This  result  is  trivial 
because  the  e-pinch  is  stable  for  arbitrary  pressure-profiles 
including  those  for  which  6=1.   As  far  as  ideal  MHD  sta- 
bility is  concerned,  this  would  then  be  the  very  comforting 
end  of  the  investigation,  but  it  is  well-known  that  life  is 
not  that  simple  and  the  unavoidable  end-losses  present  in  a 
linear  e-pinch  confirm  this  conventional  wisdom.   In  order 
to  eliminate  these  losses  one  bends  the  e-pinch  into  a  torus, 
discovers  that  there  is  no  toroidal  equilibrium  with  a  longi- 
tudinal field  alone,  hence,  introduces  a  poloidal  field  com- 
ponent, and  finally  discovers  that  a  vast  area  of  new  possi- 
bilities has  been  opened  up,  but  that  the  nice  stability 
properties  of  the  linear  e-pinch  have  been  lost  completely. 
The  reason  is  that  the  poloidal   field  necessary  for  toroidal 
equilibrium  is  so  large  that  it  introduces  rapidly  growing 
m  =  1  instabilities  for  most  of  the  configurations. 

It  is  clear  that  a  stability  analysis  of  a  diffuse  linear 
pinch  which  is  undertaken  with  an  eye  on  the  application  to 
toroidal  devices  has  to  take  into  account  the  most  obvious 
toroidal  effects  of  equilibrium  and  stability.   These  effects 
are  twofold.   Firstly,  the  condition  for  toroidal  equilibrium 
puts  a  lower  bound  on  the  magnitude  of  B  .   Secondly,  the 
finite  length  of  the  machine  introduces  a  discretization  of 


the  possible  wavenumbers  k  in  longitudinal  direction, 
because  the  wavelength  has  to  fit  into  the  torus.   The  first 
of  these  is  unfavourable  for  stability,  because  it  involves  a 
departure  from  the  stable  e-pinch  profiles,  whereas  the  latter 
is  favourable  because  it  eliminates  a  number  of  possible 
unstable  modes.   The  question  of  stability  thus  becomes  a  del- 
icate  optimization  problem  where  the  profiles  are  the 
variables  and  toroidal  equilibrium  conditions  act  as  constraints. 

On  the  basis  of  these  considerations  the  classification  of 
diffuse  linear  pinch  profiles  given  in  Fig.  1  can  be  sub- 
divided into  two  classes,  viz.  those  for  which  B   <<  B  where 
toroidal  effects  are  important  (A,B,C) ,  and  those  for  which 
Bg  _  B  where  toroidal  effects  may  be  neglected  (D,E).  Strictly 
speaking,  the  stability  theory  of  the  diffuse  linear  pinch  is 
only  a  good  approximation  for  the  latter  class  of  profiles, 
whereas  the  first  class  requires  a  consistent  toroidal  study 
of  both  equilibrium  and  stability.   Such  a  study  is  not  under- 
taken here  and  has  not  been  undertaken  yet,  to  our  knowledge. 
Instead,  we  just  introduce  the  mentioned  most  obvious  effects 
of  toroidal  equilibrium  and  stability  as  extraneous  constraints 
in  the  stability  theory  of  the  diffuse  linear  pinch.   There  is 
no  doubt   that  the  results  obtained  this  way  will  undergo 
substantial  modification  in  the  correct  consistent  treatment, 
but  we  believe  that  an  understanding  of  the  complex  influence 
of  the  profiles  on  the  stability,  already  present  in  the  simple 
straight  geometry,  is  imperative  to  progress  for  more  complicated 
geometries. 
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The  limitation  on  B  posed  by  the  toroidal  equilibrium 
requirements  is  taken  from  the  literature  by  the  method  of 
averaging  the  results  of  different  authors.   There  seems  to 
be  agreement  ^'       ''       'on   the  fact  that  fundamental  limita- 
tions on  6  as  far  as  toroidal  equilibrium  is  concerned  do  not 
exist,  but  that  a  number  of  practical  considerations  pose  a 
limit  which  is  conveniently  written  as 

Bq  <  A,  (10) 

where  A  is  the  aspect  ratio  of  the  torus,  A  =  R/a,  and 

6  max 
This  result  seems  plausible  enough  in  as  far  as  it  expresses 
that  the  field  B  ,  which  is  compressed  against  the  wall, 
should  be  large  enough  to  counterbalance   the  expansion  force 
of  the  plasma.   It  is  clear,  however,  that  it  is  no  more  than 
an  order  of  magnitude  estimate  of  an  effect  which  is  just  as 
profile-dependent  as  the  stability  results  we  are  reporting. 
We  simply  take  this  limit  as  given,  but  notice  that  the  stability 
results  improve  to  the  extent  that  one  is  able  to  keep  the 
plasma  in  equilibrium  with  a  smaller  poloidal  field  than  that 
given  by  Eq.  (10) . 

The  main  stabilizing  effect  of  the  toroidal  geometry  is 
the  possibility  to  eliminate  m  =  1  modes  with  wavelengths  larger 
than  the  circumference  of  the  torus.   This  effect  is  known  under 
the  name  Kruskal-Shafranov  limit  and  requires  a  short  exposition. 
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model  of  a  homogeneous  plasma  separated  from  a  vacuum  region 
by  surface  currents  creating  a  y-profile  as  shown  in  Fig.  3. 
For  6  <<  1  the  m  =  1  modes  have  a  growth  rate  given  by 

.   'A'  (^2) 

so  that  instability  occurs  for  |k|  <  | y (r  ) | .   This  mode  is 
eliminated  if  the  first  value  of  k  compatible  with  toroidal 
periodicity  falls  outside  the  unstable  region,  i.e.  if 
k  =  1/R  >  I p (r  ) I  or  I y (r  )r|  <  1.   In  terms  of  the  variable 
q  =  (uR)"-"-: 

q(r^)  >  1.  (13) 

The  simplicity  of  this  result  suggests  a  general  validity, 
which  can  be  rejected,  however,  by  means  of  the  same  model. 
For  6=1  the  growth  rate  is  given  by 

||£,  (14) 

so  that  instability  occurs  for  |k|  <  2 1 y (r  )  |  and  the  criterion 
for  the  absence  of  m  =  1  instabilities  becomes 

q(r^)  >  2.  (15) 

Hence, the  idea  behind  the  derivation  is  general,  but  the  explicit 
form  of  the  criterion  itself  is  not. 

The  sharp-boundary  model  with  a  vacuum  outside  the  plasma 
has  some  additional  peculiarities  illustrating   that  results 
obtained  from  special  models  can  hamper  an  understanding  of  the 
general  features.   Take,  e.g.,  a  low  6  sharp  boundary  plasma 
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surrounded  by  an  ideally  conducting  wall.   One  shows  that 
for  reasonable  wall  distances  the  plasma  is  stable  with 
respect  to  all  modes  except  m  =  1.   Observing  the  Kruskal- 
Shafranov  limit  (13)  would  then  lead  to  a  completely  stable 
configuration  and  one  can  even  show  that  fairly  high  values 
of  6  are  permitted.   However,  if  one  now  tries  to  depart 
somewhat  from  the  sharp  boundary  model,  Fig.  3  clearly  shows 
that  one  ends  up  with  a  u-profile  having  a  maximum  at  a  posi- 
tion where  the  pressure  gradient  is  large.   Suydam's  criterion 
is  violated  for  such  a  configuration  and  the  plasma  becomes 
unstable  with  respect  to  rapidly  growing  m  =  2,3,...  .  insta- 
bilities.  For  this  reason,  our  classification  of  o-stable 
profiles  of  Sec.  4  will  not  list  any  configuration  close  to 
this  model. 

VJe  now  move  to  another  class  of  special  models,  viz. 
Tokamak  configurations  (Fig.  1,A) ,  and  encounter  another 
version  of  the  Kruskal-Shafranov  limit: 

l^lmin  >    ^      °-   l^lmax^  ^  1-  <^^^ 

This  criterion  is  equally  not  general.   That  it  is  not  neces- 
sary for  stability  of  m  =  1  modes  is  proved  by  the  existence 
of  completely  stable  B  -reversed  pinch  configurations 

(Fig.  1,D),  which  have  infinite  values  of  Iplm^x'   "^^^^  ^^ 
is  not  sufficient  is  shown  by  the  fact  that  a  constant- 
pitch  region  around  the  axis  is  subject  to  m  =  1  instabilities 
having  k-values  in  a  region  around  k  =  -y (o) ,  so  that  unstable 
modes  with  |k|  >  \v\  are  not  eliminated  by  the  condition 

(16).   The  first  example  involves  an  increasing  y-profile, 
whereas  the  second  one  violates  Suydam's  criterion  on  axis. 
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Maybe  then,  the  criterion  is  a  general  statement  about 
configurations  having  a  monotonically  decreasing  y-profile 
and  satisfying  Suydam's  criterion.   There  is  a  counter- 
example for  this  statement  also,  which  is  shown  in  Fig.  4. 
Depicted  is  the  p-profile  for  a  B  -reversed  pinch  where 
p  is  a  decreasing  function  of  r,  passes  through  zero,  and 
then  takes  on  very  large  negative  values.   The  a-unstable 
region  of  k-values,  corresponding  to  radial  positions  where 
k  =  -u(r),  is  also  indicated  and  is  only  slightly  smaller 
than  the  0-unstable  region  of  k-values  which  extends  up  to 
k  =  -p(o)  on  one  side,  but  certainly  not  up  to 


the  unstable  m  =  1  modes  by  toroidal  periodicity.   It  might 
thus  be  that  criterion  (16)  is  the  proper  generalization  of 
the  Kruskal-Shafranov  limit  for  configurations  with  a  mono- 
tonically decreasing  y-profile,  not  crossing  the  axis,  and 
satisfying  Suydam's  criterion  (i.e.,  solely  for  Tokamak-like 
profiles  indicated  in  Fig.  1,  A).   It  is  certainly  a  necessary 
criterion  for  this  restricted  class  of  configurations,  because 
all  m  =  1  modes  with  |k|  <  |y(0)  |  are  unstable,  as  has  been 
proved  by  Shafranov  '   and  Robinson  ^'.      However,   we  have 
not  been  able  to  exclude  the  possibility  that  m  =  1  instabilities 
also  might  occur  for  |k|  >  |y(0)  |,  which  would  spoil  the  suffi- 
ciency of  criterion  (16) . 

For  our  purpose  the  only  useful  form  of  the  toroidal 
periodic  condition  is  the  following,  completely  general  but 
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Implicit,  condition  for  the  absence  of  the  m=l  instabilities. 
One  easily  shows  that  the  pinch  is  stable  for  fixed  m  as  k  ->  oo. 
Consequently,  the  interval  of  unstable  k-values  has  to  be 
limited: 

k^  <  k  <  k^  .  (17) 

Writing  k  =  Max  [Ik^  | ,  |k-j_  |  ^,  the  proper  generalization  of  the 
Kruskal-Shafranov  limit  then  becomes 

1/R  >  k,   or  A  <  (ka)"-^  .  (18) 

This  result  applies  to  all  configurations  and,  most  important, 
is  also  valid  in  a  a-stability  analysis.   The  interval  (k^,  k-j^) 
may  contain  stable  subintervals  (and  we  have  examples  where 
this  is  so),  but  the  possibility  of  exploiting  those  is  not 
considered  here.   However,  they  give  rise  to  so-called  stability 
windows  which  play  a  significant  role  for  the  low- shear  case. 
The  fact  that  the  criterion  is  implicit  is  no  drawback  for  the 
numerical  analysis,  because  one  has  to  check  for  the  stability 
of  m = 1,  all  k  anyway.   Of  course,  the  width  of  the  unstable 
region  (k  ,k,)  is  only  of  the  order  |m,|. 

In  conclusion,  the  consideration  of  simple  toroidal  equi- 
librium and  stability  requirements  expressed  by  the  inequali- 
ties (10)  and  (18)  leads  to  two  entirely  different  classes  of 
a-stable  configurations.   Either  m=l  is  unstable  and  one  has 
to  optimize  the  configurations  such  that  p  is  as  large  as 
possible,  but  that  still  the  condition 


Pg  <  A  <  (ka)"^  (19) 

is  satisfied.   Or  m  = 1  is  stable  and  toroidal  effects  do  not 
present  a  serious  restriction  on  the  possible  configurations. 
These  two  possibilities  correspond  to  the  cases  (A, B, C)  and 
(D, E)  indicated  in  Fig.  1. 
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4.    Classification  of  a-stable  configurations. 

The  discussion  of  the  preceding  section   leads  to  the 
following  classification: 

-  m  =  1  unstable  configurations,  stabilized  by  toroidal 
periodicity,  B   <<  B  : 

A.  Tokamak  (Figs.  lA,5-8),  decreasing  y, 

B.  B  -reversed  pinch  (Figs.  IB, 9, 10),  decreasing  u  with 
reversal, 

C.  Screw-pinch  (Figs.  IC, 11-12),  increasing  v   flattening 
in  outer  region; 

-  in=lstable  configurations,  stabilized  by  high  shear,  B  ^  B  : 

D.  Combined  pinch  (Figs.  1D,13),  increasing  y, 

E.  B  -reversed  pinch  (Figs.  lE,14),  increasing  y  with 
reversal. 

Each  of  these  configurations  will  be  discussed  on  the  basis 
of  the  numerical  results. 

The  Figs.  5-14  show  a  selection  of  o-stable  configurations, 
where  we  have  strived  to  maximize  6,  observing  the  restrictions 
posed  by  toroidal  equilibrium.   For  each  configuration  we  have 
shown:   a)  the  basic  characterizing  profile  y(r),  normalized 
for  a  wall  radius  a  =  1,  b)  the  functions  p (r) ,  Bg  (r) ,and 
B  (r) , which  determine  the  equilibrium  completely,  c)  the 
derived  current  profiles  JQ(r)  and  i^i^)    which  are  of  practical 
interest.   We  indicate  the  values  of  6  and  B„  as  defined  by 
the  Eqs.  (2)  and  (11) ,  and  the  value  of  ^  ^  as  defined  by 
Eq.  (8).   For  the  configurations   that  are  m  =  1  unstable  we 


also  give  the  maximum  growth  rate  ~  ^  =x  °^  ^^^   m  =  1  insta- 
bilities and  the  maximum  aspect  ratio  A    for  which  criterion 

max 

(18)  is  still  satisfied. 
A.    Tokamak. 

The  sequence  is  started  with  two  configurations  (Figs. 5,6) 
which  exceed  the  limit  (10)  posed  by  toroidal  equilibrium  by 
some  orders  of  magnitude.   They  have  high  values  of  B  and  the 
poloidal  current  is  larger  than  the  toroidal  current.   The 
purpose  of  these  pictures  is  to  show  the  smooth  transition 
from  9-pinch  to  o-stable  Tokamak  profiles.   If  one  could 
keep  these  configurations  in  equilibrium  by  means  other  than 
increasing  the  poloidal  field,  one  would  have  a  high- B  stable 
Tokamak  with  the  additional  advantage  of  enhanced  heating  by 

the  poloidal  current.   These  advantages,  combined  with  the 

12  ) 
finding  of  Callen  and  Dory     "that  from  magnetohydrodynamic 

equilibrium  considerations  alone  there  is  no  apparent  limit 

2 
on  B   up  to  about  A   and  thus  limit  on  B  up  to  about  unity  , 

warrant  a  further  study  of  the  practical  limitations  on 
toroidal  Tokamak  equilibrium  listed  in  the  mentioned  reference. 

The  next  two  pictures  (Figs. 7, 8)  show  more  standard 
Tokamak  profiles.   These  represent  a-stable  configurations 
with  the  highest  values  of  6  we  could  obtain  without  violating 
condition  (10)  for  toroidal  equilibrium.   The  profiles  are 
the  result  of  a  subtle  optimization  process,  where  the  condi- 
tion (10)  for  toroidal  equilibrium  favors  a  large  value  of 
u  at  the  wall,  the  toroidal  periodicity  condition  (18)  for 


the  absence  of  m  =  1  modes  favors  a  small  value  of  y  on  axis, 
but  the  stability  conditions  for  a  diffuse  linear  pinch  in 
general  demand  a  large  shear,  i.e.,  the  difference  u(o)  -  U (a) 
as  large  as  possible.   For  a  Tokamak,  where  the  p-profile 
is  decreasing  without  crossing  the  axis,  these  conflicting 
conditions  lead  to  severe  limitations  of  the  possible  values 
of  6  and  one  ends  up  with  values  of  6  of  the  order  0.03  to 

0.04.   The  demonstration  that  Tokamak  profiles  exist  which 

l4) 
are  stable  with  respect  to  m  =  2,3,...  is  due  to  Shafranov 

In  order  to  satisfy  Suydam's  condition  at  r  =  0   the 

pressure  profile  has  to  be  rather  flat  on  axis,  such  as  is 

4 
given  by  p ~ exp (-Br  ).   A  consequence  of  this  condition  and 

the  equilibrium  equation  (1)  is  the  structure  of  the  current 

profile  jg,   which  tends  to  develop  a  reversal  of  the  sign. 

More  smooth  j   profiles  can  be  obtained  by  starting  from  a 

2    4 
pressure  profile  like  p  •^  exp (-Ar  -Br  ),  but  the  result  is 

evidently  Suydam-unstable  plasma  on  axis  (See  Sec. 5,  Fig. 15). 
B.    B  -reversed  pinch. 

In  view  of  the  mentioned  limitations  on  the  permissible 
6  in  a  Tokamak,  originating  as  a  result  of  conflicting  condi- 
tions on  the  vi-profile,  a  logical  step  seems  to  consider  mono- 
tonically  decreasing  y-profiles  which  reverse  sign.   The  first 
configuration  (Fig. 9)  shows  that  this  can  be  done  by  inducing 
currents  of  opposite  direction  through  the  outer  regions  of 
the  pinch.   The  result  is  a  y-profile  with  very  favourable 
properties  for  magnetohydrodynamic  stability,  viz.  rather  small 


values  of  y  on  the  average  and  still  high  shear.   Accordingly, 
one  can  obtain  stable  configurations  with  6  =  0.3,  which  is 
an  order  of  magnitude  larger  than  the  possible  values  for  a 
Tokamak.   The  m  =  1  modes  are  still  unstable,  because  v    is 
a  decreasing  function  (see  Sec. 3),  but  they  are  eliminated 
by  means  of  the  toroidal  periodicity  condition,  whereas  the 
higher-m  modes  are  stabilized  by  the  presence  of  high  shear. 
Toroidal  equilibrium  is  due  to  the  large  poloidal  field  at 
the  wall.   The  configuration  can  be  considered  as  a  super- 
position of  Tokamak  and  screw  pinch. 

From  a  practical  point  of  view  the  presence  of  large 
currents  close  to  the  wall  seems  unwanted.   One  can  cure 
this  by  continuing  the  current  profiles  such  that  they  decay 
to  zero.   In  case  the  B  -field  is  partly  created  by  external 
currents,  this  will  lead  to  a  minimum  in  the  p-profile  because 
the  absolute  value  of  y  for  a  vacuum,  field  has  to  decay  like 

/r  .   The  associated  instabilities  have  small  growth  rates 
so  that  a-stability  is  still  possible,  but  the  value  of  6  has 
to  be  reduced  (see  the  similar  effect  for  the  screw  pinch  in 
Fig.  12).   Another  possibility  is  shown  in  Fig.  10.   Here, 
there  is  no  minimum  in  y  because  the  B  -field  is  created 
exclusively  by  means  of  plasma  currents,  so  that  B   -*-  0  and 
y  ->  "  at  the  wall.   The  resulting  high  shear  is  favourable 
both  for  stability  and  toroidal  equilibrium,  so  that  the  large 
values  of  6  can  be  maintained.   The  configuration  is  m  =  1 
unstable,  but  we  have  already  shown  in  the  preceding  section 
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that  the  unstable  region  occurs  for  k-values  which  can  be 
eliminated  by  means  of  the  toroidal  periodicity  condition 
despite  of  the  large  value  of  u  at  the  wall. 

The  configuration  of  Fig.  10  can  be  looked  upon  as  a 
superposition  of  Tokamak  and  combined  pinch.   The  sequence 
of  a-stable  profiles  could  be  continued  further  with  super- 
positions of  Tokamak  and  B  -reversed  pinch,  having  then  a 
u-profile  going  through  zero  and  infinity,  but  the  increas- 
ing complexity  of  the  current  distributions  (and,  therefore, 
the  associated  necessity  of  complex  circuit  programming) 
makes  such  configurations  of  little  practical  use.   From 
the  completely  different  distributions  of  Figs.  9  and  10, 
it  is  clear  that  many  more  stable  configurations  with  high 
values  of  B  can  be  constructed.   To  our  knowledge,  stable 
Bg-reversed  pinches  have  not  been  considered  before,  neither 
theoretically  nor  experimentally. 
C.    Screw-pinch. 

Another  sequence  of  a-stable  profiles  is  obtained  for 
increasing  y-profiles.   Such  configurations  arise  as  a  result 
of  current  distributions  which  are  completely  different 
from  those  of  a  Tokamak.   In  a  Tokamak,  the  toroidal  current 

has  to  be  fairly  peaked  on  axis  in  order  to  create  a  y-profile 

l4 ) 
with  enough  shear  for  stability  ' .      In  screw  pinches,  com- 
bined pinches,  and  B  -reversed  pinches  the  necessary  shear 
is  created  in  exactly  the  opposite  manner,  viz,  by  current 
densities  that  are  small  on  axis  and  large  in  outward  regions. 


24 


Such  distributions  arise  as  the  result  of  imperfect  induction 
of  skin-currents  (screw-pinch)  or  deliberate  programming 
(B  -reversed  pinch) . 

The  stability  of  screw-pinch  configurations  was  considered 
first     for  a  sharp-boundary  model  with  a  step-function  for 
y (r) .   Unlike  the  sharp-boundary  model  with  a  vacuum  outside 
the  plasma  column  (see  Sec. 3),  this  model  can  be  turned  into 
an  acceptable  diffuse  configuration  where  the  u-profile  is 
monotonically  increasing.   Such  a  profile  is  more  or  less 
the  mirror-image  of  that  of  a  Tokamak  (cf.  Figs. 7  and  11), 
so  that  a  similar  optimization  game  can  be  played  as  before, 
but  with  slightly  more  favourable  rules.   The  stability  re- 
quirements on  axis  are  now  automatically  satisfied  because 
U  is  small  there.   Toroidal  equilibrium  and  linear  pinch 
stability  require  a  large  value  of  v    at  the  wall,  whereas 
the  toroidal  periodicity  condition  for  the  absence  of  m  =  1 
instabilities  favors  small  values  of  y (a) .   These  conflicting 
conditions  again  restrict  the  possible  values  of  B,  but  not 
to  the  same  extent  as  for  a  Tokamak,   Fig.  11  shows  that 
stable  configurations  with  values  of  6  up  to  0.3  can  be 
obtained. 

The  unattractive  feature  of  large  current  densities  at 
the  wall  can  be  remedied  as  shown  in  Fig.  12.   Because  the 
fields  have  to  transform  into  a  vacuum  configuration  the 
function  u(r)  now  gets  a  maximum,  so  that  Suydam's  criterion 
is  no  longer  satisfied.   Since  the  pressure  gradient  is  also 
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small  in  that  region  the  growth  rates  of  the  instabilities 
are  still  small,  so  that  a  a-stable  configuration  is  still 
possible.   However,  the  value  of  6  has  to  be  reduced  to 
about  0.08. 

The  necessity  to  observe  the  toroidal  periodicity  con- 
dition (18)  for  these  configurations  arises  as  a  result  of 
the  fact  that  the  shear,  especially  in  the  outer  region  of 
the  pinch,  is  insufficient  to  stabilize  the  m  =  1  mode. 
This  is  even  more  salient  for  the  case  that  the  y-profile 
has  a  maximum  like  in  Fig.  12.   For  this  case,  a  proof  has 
been  given    for  the  instability  of  the  m  =  1  mode,  even  if 
p'  =  0  (so,  it  is  a  distinct  phenomenon  from  violation  of 
Suydam's  criterion),  similar  to  the  proof  for  monotonically 
decreasing  p-profiles.   The  growth  rates  of  these  instabili- 
ties are  quite  large  (cf.  the  values  of   -fi   in  Figs.  11 
and  12  with  table  I) ,  so  that  it  is  imperative  to  satisfy 
condition  (18) .   An  alternative  cure  of  these  instabilities 
is  to  enhance  the  shear  by  increasing  B  and  decreasing  B  in 
the  outer  regions  of  the  pinch.   This  leads  to  the  next  class 
of  profiles. 
D.    Combined  pinch. 

A  completely  new  class  of  configurations  is  obtained 
if  one  increases  the  shear  and  the  ratio  Bq/B  such  that 
the  m  =  1  modes  become  stable.   For  these  configurations 
toroidal  equilibrium  and  stability  considerations  can  be 
neglected  and  one  has  to  deal  with  the  stability  of  the 
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diffuse  linear  pinch  alone.   A  continuous  transition  from 

screw-pinch  to  combined  pinch  and  B  -reversed  pinch  has  been 

described  in  Ref.  4  by  means  of  a  sharp  boundary  model  with 

force-free  currents  outside  the  plasma  column.   If  one 

insists  upon  having  vanishing  current  densities  at  the  wall, 

this  transition  is  much  less  continuous  and  one  finds  that 

only  one  type  of  configuration  survives  which  has  a  monotoni- 

cally  increasing  u-profile  is  stable  for  m  =  1,  and  has  no 

current  at  the  wall.   This  is  the  combined  pinch^^ '  where 

the  B  -field  tends  to  zero  at  the  wall  so  that  v  (a)  ->■   «>. 
z 

The  resulting  configuration  shown  in  Fig.  13  is  one  of  the 
nicest  we  encountered.   It  is  a  pure  pinch  configuration 
where  confinement  is  solely  due  to  the  interaction  of  plasma 
currents  and  the  wall.   In  the  equilibrium  state,  no  external 
current  is  present.   The  very  large  shear  stabilizes  all  modes, 
including  m  =  1.   The  current  distribution  is  naturally  ob- 
tained by  smoothing  a  sharp  boundary  stabilized  z-pinch  with 
a  trapped  B  -field,  surface  currents  and  no  B  -field  outside 
the  plasma  column  (i.e.,  a  step  function  0  -  <»  for  u  (r)  )  . 
On  the  basis  of  the  latter  model,  the  possiblity  of  stable 
confinement  of  plasma  with  high  values  of  6  has  been  realized 
as  early  as  1957-'-®^ 
E.    B  -reversed  pinch. 

The  large  shear  of  the  combined  pinch  can  even  be  ex- 
ceeded if  one  allows  for  negative  values  of  the  B  -field  in 
the  outer  region  of  the  pinch.   The  resulting  B  -reversed 
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configurations  ''      ^''      ^',    shown  in  Fig.  l4,  has  the  highest 
possible  shear  in  the  sequence  of  configurations  with  mono- 
tonically  increasing  ji-profiles,  like  the  Bg-reversed  pinch  in 
the  sequence  of  raonotonically  decreasing  M,-profiles.   The 
possible  values  of  P  are,  therefore,  comparable  for  both  field- 
reversed  pinches,  i.e.  0.2  -  0.3-   Like  in  the  combined  pinch, 
the  current  profiles  nicely  decay  to  zero  at  the  wall,  so  that 
the  outermost  region  has  a  vacuum  configuration.   In  the 
equilibrium  state,  the  reversed  B  -field  is  maintained  by  an 
external  poloidal  current  of  opposite  sign  with  respect  to  the 
plasma  jg-current.   As  shown  in  Ref.  15,  the  wall  position  and 
the  amount  of  field  reversal  determine  the  stability  in  such  a 
way  that  the  allowable  percentage  of  the  field  reversal 
decreases  if  the  wall  position  is  increased,  and  vice  versa. 
Therefore,  the  sequence  of  configurations  with  decreasing  ii- 
profile  naturally  ends  here,  because  the  introduction  of  more 
shear  apparently  is  not  possible  without  at  the  same  time 
introducing  rapidly  growing  current-driven  instabilities. 

The  given  profiles  in  Fig.  l4  differ  in  one  important 
aspect  from  the  ones  studied  in  Ref.  15,  which  have  a  hollow 
pressure  profile  and  a  j  -current  profile  which  is  peaked  on 
axis.   Here,  the  features  are  Just  reversed,  viz.  the  pressure 
has  a  maximum  on  axis  and  the  j  -current  has  a  maximum  off- axis. 
Although  both  configurations  have  comparable  stability  proper- 
ties, one  of  them  could  be  a  preferred  equilibrium  arising  in 

20) 
the  process  of  creation  and  diffusion  of  the  plasma.   ' 


5.    Unstable  configurations. 

In  order  to  arrive  at  a  balanced  conclusion,  it  seems 
appropriate  also  to  expose  the  other  side  of  the  picture. 
For  that  purpose,  we  have  chosen  three  representative  unstable 
configurations  which  can  well  be  distinguished  from  their 
stable  counterparts,  but  are  close  enough  to  feel  uncomfortable. 

In  Flgo  15  an  unstable  Tokaunak  configuration  is  shown 
together  with  a  stable  one.   The  unstable  configuration  deviates 
from  the  stable  one  in  that  the  pressure  decays  faster  on  axis 
aind  the  pololdal  field  is  more  concentrated  towards  the  axis, 
but  the  deviation  does  not  seem  to  be  very  large  until  one 
considers  the  more  detailed  information  given  in  the  derived 
quantities  [i,    j'  ,  j   and  the  Suydam  function  S  =  1  +  -^^  {^) 

given  in  Fig.  I6.   It  turns  out  that  the  instability  is  caused 
by  violation  of  Suydam' s  criterion  on  axis,  so  that  there  is 
Instability  for  all  m  (except  m=l  which  is  still  absent  because 
of  toroidal  periodicity).   The  growth  rates  are  very  large. 

In  Fig.  17  an  unstable  screw-pinch  is  shown  together  with 
a  stable  one.   Again,  the  profiles  are  rather  close,  but  the 
derived  quantities  (Fig.  18)  give  much  more  information.   For 
a  screw-pinch  configuration  violation  of  Suydam' s  criterion  on 
axis  would  not  matter  very  much  because  \i   is  small  there,  so 
that  the  growth  rate  is  small  and  a  large  number  of  higher  -m 
modes  are  eliminated  by  toroidal  periodicity.   However, 
instability  is  now  induced  by  decreasing  the  current  density 
close  to  the  wall  so  that  a  region  of  small  shear  develops. 
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Hence,  although  Suydam's  criterion  is  satisfied  throughout  the 
plasma,  het  higher-m  modes  are  unstable  having  rather  large 
growth  rates o 

Finally,  we  present  in  Fig.  19  an  example  which  specifi- 
cally has  been  constructed  for  the  purpose  of  showing  the 
intricacies  of  the  problem.   We  show  an  m  =  1  unstable  combined 
pinch  configuration  which  is  extremely  close  to  a  a-stable  one 
but  still  has  very  large  growth  rates.   From  the  equilibrium 
profiles  in  Fig.  19  it  would  even  be  impossible  to  predict  which 
of  the  two  is  the  stable  and  which  is  the  unstable  configura- 
tion.  Figure  20  reveals  how  this  example  was  constructed.   The 
current  profiles  were  provided  with  a  small  dent  such  that  a 
tiny  plateau  in  the  |j,-profile  was  produced.   Consequently,  in 
that  small  region  Suydam's  criterion  is  violated  and  the  con- 
figuration becomes  unstable  with  respect  to  a  gross  m  =  1  mode 
having  some  tendency  to  localize  in  the  Suydam-unstable  region. 
This  shows  that  a-stable  bands  of  equilibrium  profiles  only 
occur  for  profiles  that  are  close  with  respect  not  only  to  p, 
Bg,  and  B  ,  but  also  with  respect  to  jg  and  j  . 

Whether  dissipative  effects  provide  sufficient  smoothing 
of  the  profiles  in  order  to  prevent  the  occurrence  of  situa- 
tions like  just  described  is  an  open  question.   On  the  other 
hand,  the  eigenfunctlons  also  might  be  required  to  be  suffi- 
ciently smooth.   Numerically,  artificial  smoothing  of  the 
eigenfunctlons  can  easily  be  obtained  by  putting  an  upper  bound 
on  the  allowable  values  of  the  radial  derivative  of  ^  ,  thus 
eliminating  local  instabilities.   Such  a  procedure  enlarges 
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the  class  of  o-stable  profiles  significantly  and,  moreover, 
produces  broad  bands  of  a-stable  profiles.   This  heuristic 
procedure  might  find  a  justification  in  a  more  general  descrip- 
tion of  the  plasma,  incorporating  f inite-Larmor  radius  effects. 
The  given  examples  suffice  to  show  that  categoric  answers 
about  which  configurations  are  bad  and  which  are  good  with 
respect  to  stability  cannot  be  given. 
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6.    Conclusion. 

On  the  basis  of  the  a- stability  theory  of  the  diffuse 
linear  pinch,  supplemented  with  simple  toroidal  equilibrium 
and  stability  requirements,  a  classification  of  a-stable  high- 
shear  configurations  has  been  given,  which  we  believe  is  fairly 
complete o   Five  different  groups  of  configurations  are  obtained 
characterized  by  the  function  |J.(r).   In  each  of  these  groups  a 
number  of  a-stable  configurations  exist  which  have  rather  high 
values  of  p.   However,  highly  unstable  configurations  with 
equilibrium  profiles  very  similar  to  the  stable  ones  also  occur 
(Sec.  5)  and  it  seems  very  difficult  to  avoid  these  unstable 
states.   It  is  possible  to  refine  the  description  of  the  plasma 
considerably,  especially  by  taking  toroidal  curvature  and 
resistive  effects  into  account,  but  the  mentioned  problem  will 
remain.   Its  solution  will  require  a  dynamical  theory  of  the 
pinch  abolishing  the  artificial  separation  of  plasma  equilibrium 
and  stability. 
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Fig.  1.   Classification  of  linear  pinch  configurations  accord- 
ing to  u-profile.   A.   Tokamak;   B.   B^-reversed  pinch; 
C.   Screw-pinch;   D.   Combined  pinch;   E.   B^-reversed 
pinch. 
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Fig.  2.   Striae  of  a-stable  profiles,  3  =  0.10,  Sg  =  0.23  and 
^  2  ^  4  X  IQ-^^ 
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Fig.  3.   Sharp-boundary  model  for  a  pinch  surrounded  by  vacuum. 
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^  o 


Fig.  4.   y-profile  for  a  B  -reversed  pinch.   In  box  is  the 
detail  of  y-profile  near  origin. 
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Fig.  5.     Non-equilibrium  Tokamak,  6  =  0.67, 
2 


=  2  X  10  , 


m  =  1  unstable. 
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Fig.  6.   Non-equilibrium  Tokamak, 


.  0.20,  6,  •  480,  A„^-  5.0, 
9  X  10~  ,  m  =  1  unstable. 
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Fig.    7.      Toka:aak,    6   =    0.037,    gg    =    1.73,    A^^^   -    2.17, 


=   4   X   lo'-'-^,    n^   =   -    1.7   X   lO'^,    m  =   1   unst^^ble. 


.60n 


Fig. 


Tokamak,    B    =    0.037,    6      =    1.43,    A  =   1.89,    JZ  ''^10 
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Fig.  9.   Bg-reversed  pinch,  B  =  0.30,  gg  =  2.28,  A^^  =  2.52, 
^  •^  =  4  X  10~^^,  fi^^   =  -  lO"'',  m  =  1  unstable. 
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iGh,    6   =    0.30,     6„    =    0.38,"    A  =    1.40, 


at   the  wall,    in  =   1   unstable. 


45 


Fig. 11.      Screv7-pinch,    B   =   0.30,    B^   =   1.71,    A^^,^  =   1.89, 

)~"  "    =    4   X   lO""*"^,    i'?' ,      =   -   7.0  X   lO"^,    IP.  ■■=   1  unstable 
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Fig. 12.   Screw-pinch,  P  =  0.08,  Pg  =  2.70,  \^   =  ^-50, 

T— 2  _  ,  -.  10"'''  n^   =  -3.0  X  10"  ,  no  current  at  the 
wall,  m  =  1  unstable. 
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Fig. 13.      Combined  pinch,    p  =   0.10, 
m  =   1   stable. 


0.22,    T~     =  4   X   10" 


Fig. 14.  B  -reversed  pinch,  P  =  0.20,  Pg  =  0.38,  YZ  =  '^  ^  10 
m  =  1  stable.  It  becomes  unstable  for  wall  position 
farther  than  0.80a. 
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Fig. 15,   Unstable  Tokamak  (dotted  lines)  compared  with  stable 
one  (solid  lines),   p  =  O.036,  Pg  =  I.30,  stable  one: 
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Fig. 16.   Profiles  n,  jg,  j^,  and  Suydam  function  S  for  Tokamak 
configurations  of  Fig.  15' 
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Fig. 17.   Unstable  screw-pinch  (dotted  lines)  compared  with  stable 
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Fig. 18.   Profiles  u,  5q.    3^^    ^&   S  for  screw-pinch  configurations 
of  Fig.  17. 
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Fig. 19.   Unstable  combined  pinch  (dotted  lines)  compared  with 


stable  one  (solid  line),   p  =  0.10, 


0.23,    stable 


one:      5^       =  4   x   lo"      ,    m  =   1   stable,    unstable    one: 


3   X  10' 


=   1). 
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Fig. 20.   Profiles  n,  Jg,  j^  and  S  for  combined  pinch  configurations 
of  Fig.  19. 
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